Temperature effects on quantum cloning of states and entanglement 
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Introduction. — It has been known that based on the prin- 
ciples of quantum mechanics, accurate cloning of arbitrary 
quantum state is impossible JU]. Nevertheless, on the account 
of the significant role of copying in quantum computation and 
quantum communication, a variety of approximate quantum 
cloning transformations have been proposed, e.g., universal 
cloning (UC) machines producing identical copies from ar- 
bitrary input qudits 0, B 0, IH la], phase-covariant cloning 
(PCC) machines of a class of partially known qud its 0, 0], 
and optimal entanglement cloning machines J^, I10H. For re- 
cent reviews see Refs. fiUl . 

A question of practical relevance is how uncontrollable en- 
vironmentally induced decoherence or dissipation can affect 
performance of quantum cloning machines. In closed sys- 
tems, an initially pure state evolves to another pure state. In 
practice, however, preparation of pure states and/or keeping 
them pure are not generally easy tasks. In general, interac- 
tion with an environment degrades purity of quantum systems 
and makes their states mixed. A usual effect that a thermal 
environment can cause is thermalization (as a kind of dissipa- 
tion) UUdlt]. That is, because of interaction with the environ- 
mental degrees of freedom which are in thermal equilibrium, 
the quantum system will also be driven toward equilibrium. 
It should be noted that a generic isolated quantum many-body 
system does also relax to a state well described by the standard 
statistical-mechanical prescription lfl4ll . In this paper, our aim 
is to investigate temperature effects on the performance of the 
cloning machines. It has been known that decoherence can put 
a limitation on the number of clones that a quantum doner can 
generate lfl5ll . There is also a model in which the robustness 
of the doner increases with the number of qubits 111611 . 

Through a thermalizing process, the density matrix of the 
system q in long time will approach the Boltzmann state 
£>th = e~P H jZ, where ksP is inverse temperature (fcs is 
the Boltzmann constant), H is the Hamiltonian of the open 
quantum system, and Z = Tr(e~ 3ff ) is the partition func- 
tion. Energy and phase relaxation processes, with the time- 
scales Ti and T2, respectively, are common processes present 
when approaching an induced equilibrium state. For a more 
precise and elaborate discussion of dissipation, thermaliza- 
tion, and decoherence see Ref. lfl2ll . Some (phenomenolog- 
ical) models for the underlying dynamics of the thermaliza- 
tion have already been proposed lfl3ll . We assume that the 



time-scale in which typical correlation functions of the en- 
vironment decay, i env ., is much smaller than all other time- 
scales, i.e., i env . <C T c ,Tdiss. = rnin{Ti, T2, To}, where r c is 
the time-scale of the cloning process and To is the time-scale 
dictated by all other relaxation mechanisms. This assumption 
is important for the Markovian analysis of the dynamics of the 
thermalization lfl7l [Till . This implies that during the cloning 
process, a negligible amount of information flows from the 
environment to the system (or vice versa). Here, we also as- 
sume that r c < Tdi ss .. This extra condition allows us to ig- 
nore dynamical effects of the thermalization, hence consider 
a simple static (toy) model — explained below — to bring 
temperature into play. Despite these simplifying assumptions, 
we will argue that the result is still reliable enough to give a 
hint about how temperature effects can change performance 
of different cloning machines such as the universal doners, 
phase-covariant doners, and entanglement doners. Indeed, 
such investigation has an immediate importance in attempts 
to realize quantum cloning in systems where (due to thermal 
and perhaps other noise effects) the preparation of pure states, 
whether initial state of the system to be cloned or the quantum 
hardware, is difficult, such as in NMR systems l[l9i |20T| . For 
another study using a different approach, see Refs. l2ll I22I1 . 
For the purpose of illustration, we only consider the case of 
symmetric 1 — > 2 qubit doners. Extension to M — > N qudits 
is straightforward as well. 

Optimal universal and phase-covariant cloning transfor- 
mations. — In the universal cloning transformation, it is usu- 
ally assumed that the qubit state to be cloned is a pure state, 
|*} a = cos § 1 0) + e^sinf|l), and the blank copy ib) and 
the quantum cloning machine (also called ancillary system, 
c) are each in a known pure state, say |0) |S H. S HU. 
The symmetric cloning transformation, then, acts in this 
way: U (|*) a |0) 6 |0) c ) = \T) abc , where Tr bc (|T) abc (T|) = 
Tr ac (|T) a h c (T|). The latter condition guarantees that the fi- 
nal clones both have the same states, p° ut = g™ 1 . A mea- 
sure to quantify performance of a cloning machine is the fi- 
delity between the original and the output states, F(p, a) = 



Tr^/^72 

ap 1 / 2 J . Optimization of the fidelity over all input 

states on the Bloch sphere results in the qubit optimal univer- 
sal doner, in which F = 5/6 J3.[3t]. For orbital states, where 
6 is an a priori known constant and <j> £ [0, 2tt), a class of 
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FIG. 1: Fidelity (F) of UC vs 6 for some values of r\: (a) e = 5/11 
and (b)e = 2/3. 



FIG. 2: Variation in the fidelity of PCC with 6 for: (a) e = 5/11 and 
(b) e = 2/3, and different values of rj. 



phase-co variant cloning machines has been suggested J^]. Af- 
ter the cloning process, in the computational basis {|0), |1)} 
(the eigenvectors of a z = diag(l, —1)) each of the clones can 
be identified by the density operator: £>qq = /-t 2 £?o + 1,2 an d 
0§? = 2fivg$ 1 , where /j 2 + 2v 2 = 1, and v" 1 = 1/6 for UC 

and v 1 = (1 7— i — rs)/ 4 for PCC. Most of this descrip- 

tion is also valid when the original quantum system is initially 
mixed. 

Our main assumption is that preparation of the initial pure 
state | "J) is diluted by a thermal bath in the following special 
and simple form: 

e in = (l-e)|*)(*| + e fth , 0<e<l. (1) 

The parameter e, which measures how thermally perturbed the 
preparation is, may in general be time-dependent. Nonethe- 
less, based on our earlier assumptions, it would be a fairly 
slow-varying time-dependent function so that with a good ap- 
proximation we can take it a relatively small constant of the 
order of r c /Tdi ss .. This state does not seem to arise natu- 
rally from a typical thermalization dynamics. Nevertheless, 
in Ref. 01 it has been illustrated that general behaviors ob- 
tained from such a simple preparation assumption (in the con- 
text of the geometric phases) have general features similar to 
those obtained from the Lindblad equation for the dynamics. 
It is worth mentioning that in the limit of infinite tempera- 
ture, the thermalized density matrix p t h is equivalent to pure 
noise 112411 . In that case, e represents the degree of pure noise 
existing during the process (for example, in the case of NMR 
systems, due to fluctuations of the external magnetic fields and 
similar reasons). A more general analysis of quantum cloning 
in the presence of a thermalization mechanism is yet lacking, 
but our simple analysis may also shed some light before hav- 
ing a more complete analysis at hand. 

First, we consider the effect of the thermal term only on the 
state of the doner, that is, the quantum cloning hardware is 
thermally diluted as in Eq. ([TJ. In this case, the initial state 
of the machine is mixed. Considering the fact that in the op- 
timal UC and PCC, the initial state of the cloning machine 
can be any pure state 0, HU, one can conclude the opti- 
mal fidelity here is achieved by the existing optimal cloning 
transformations. By a similar analysis, it appears that for the 
case of diluted joint blank and ancillary systems, one can con- 
sider the joint state as a new blank copy and attach some new 
reservoir to the whole Hilbert space of the input states (i.e., 
the information qubit, the blank copy, and the ancilla state) as 
a new ancillary system and then define a new transformation 



for cloning j23ll . This would in fact be the existing optimal 
cloning transformation, now acting on a larger Hilbert space, 
and hence one obtains the same optimal fidelity again. How- 
ever, from an experimental point of view, thermalization ef- 
fects are likely to occur during the cloning process rather than 
at the initial preparation level — for instance in NMR sys- 
tems lfl9ll20ll . Therefore, to be more precise, thermal effects 
during the preparation should also be taken into account. 

Now, we consider the case in which the input state |\&) is 
thermally diluted as in Eq. (Q~|i. Our aim now is to compare 
the similarity between the clones and the input state of our in- 
terest, i.e., I^). Indeed, here we assume that the model of the 
cloning machine consists of two parts: the first is the dilution 
of the input pure state which models the imperfect feature of 
the machine, and the second is some known UC or PCC trans- 
formation which is applied to this diluted state. The Hamilto- 
nian of the qubit system is taken to be H = luq(j z /2 (luq > 0), 
whence, Z = 2 cosh rj, where ?; = luq/3/2. More general 
cloning transformations in spin networks with more compli- 
cated Hamiltonians can be found, for example, in Ref. 12111 . 
The fidelity of the output state and the unperturbed initial state 
can be calculated as follows: 

F(6, e, 77) =/i 2 [l - e + e(e^ cos 2 - + e" sin 2 -)/Z] 
+ {jiv - /J 2 /2) (1 - e) sin 2 9 + v 2 . 

Figure[T]illustrates how the fidelity in the UC behaves in terms 
of 8 (orbit of the state) in thermally diluted states, for two 
different values of e (the degree of thermalization) and rj (oc 
1 /T). It can be seen that when 

9 9 
e < cosh 77 /(e~ r ' sin 2 -^+e v cos 2 -), (2) 

the fidelity of the UC is higher than the classical value 1/2. 
This threshold is the fidelity of a classical-like 1 — > M uni- 
versal cloning in which with a given probability, an unknown 
input state is sent to one of the M parties and a completely 
randomized state is transmitted to any of the other ones, of 
course, in the limit of large M l25ll . In the literature, however, 
"classical doner" has been attributed to some other cloning 
transformations as well — see l3l l26ll . In other words, in some 
cases thermal noise (even in the simple form of Eq. (Q])) can 
result in a lower performance than a classical machine. For 
9 ^ tt/2, the condition © implies that for all ^ e < 1, the 
fidelity of the output of the UC is always greater than that of 
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FIG. 3: Fidelity vs e and 9 in low temperature limit (77 — + 00): UC 
(left) and PCC (right). 



the classical doner (if u>o was negative, this would occur for 
9 ^ 7r/2). Equation (|2]i can also be interpreted as a condition 
on temperature for a given 9 and e in order to outperform a 
classical doner. Figure |2] shows the variation of the fidelity 
of the outputs of the PCC machines in terms of 9, for some 
fixed values of e and r\. As is clear from this figure, in the case 
of equatorial qubits, similar to the case of the UC, the fidelity 
of the outputs does not vary with temperature — according to 
Eq. ©, this feature is due to the symmetry property of such 
states. Low temperature limits of the fidelity for both UC and 
PCC have been depicted in Fig. [3] In the case of the UC, for 
all 9 in [0, w), the fidelity is a decreasing function of e. The 
corresponding graph for the PCC also shows a decrease in the 
fidelity for different values of 9 E [0, 7r/2) with the pertur- 
bation factor e. However, a closer inspection shows that here 
there are also some 9s (> 2.52 and less than it rad) in which 
the fidelity of the PCC is an increasing function of e. At high 
temperature limit, the fidelity of both UC and PCC, for all 9s, 
is a decreasing function of e. Another important point that can 
be concluded from the figures is that in some cases, the quality 
of the clones at the output of the UC can be better than that of 
the PCC — see for example those regions of Fig. [3]in which 
e and 9 are large and small, respectively. This is indeed con- 
trary to what happens when the cloning is performed perfectly 
without any external noise. 

Entanglement cloning. — Quantum cloning can be used to 
clone or broadcast entanglement as well [HI H3, l28ll . 
Let us assume that we have an initial state in the form of 
\^a)ab = a|01) O 6 - Vl - a 2 |10) a6 , where a is real and 
I a I ^ 1. As in the cases of the UC and the PCC, sup- 
pose that because of a thermal environment, the initializa- 
tion is diluted as in Eq. (Q}. Let us take our system to be 
two spin- 1/2 particles interacting via the XX Hamiltonian: 
H = J(er°cr* + Vy(Ty), where a x and <r y are Pauli matri- 



ces. Now, we want to compare performances of the following 
schemes of entanglement broadcasting between two parties in 
the presence of thermal noise: (i) Local cloning by means of 
two optimal UC machines copying each qubit separately l27ll . 
In this scenario, after the cloning process and discarding the 
ancillas, we will have the overall state g aa 'bb' whose two first 
(last) qubits are the copies of a (b). (ii) Non-local cloning of 
the two-qubit state as a whole with the UC machine of 4-level 
quantum states . (iii) Cloning by an optimal entanglement 
doner H. 

After some algebra, it can be seen that the density matrices 
of the clones in cases (ii) and (iii), and g a >b (also Q a y, g a b, 
and Qa'b') — nonlocal copies — in case (i), read as follows: 
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L(l-e)(2a 2 -l)]|01)(01| 
L(l-e)(2a 2 -l)]|10)(10| 



Q 0Ut 


= (^ + 


1-M> 
4 > 


1=1 + 


e cosh 7 ^ | 


1-Af 


2 ' 


Z > + 


4 


1=1 + 


e cosh 7 ^ | 


1-Af 


2 ~ 


Z > + 


4 



+[M( 
+[M( 

-M[(l - e)aVT~^ + § sinh 7 ](|01)(10| + |10) (01|), (3) 



in which 7 = 2(3 J, Z = 2(1 

\/l + 4M -9M 2 )/26, Mi = 

6A 2 + 4AC, A 



- COSI17), L = 3(1 + 2M 
(2/3) 2 , M fl - 3/5, M ffl 

A(y 



yJ{l/2 + l/VT3)/3, and C = A(Vl3 
3)/2. Note that, the output states of case (ii) for all values of 
e, a, and 7, the nonlocal copies of case (i) g a 'b, and the output 
states of case (iii) for e = 1 and V7 or a = ±1/^/2 (for all 
e and 7) all can be written in the following compact form: 
g out = Mg m + (1 - M)I/4, where / is the 4 x 4 identity 
matrix. 



To determine the regions in which the output states are 
separable or inseparable, we use the well-known Peres- 
Horodecki positive partial transposition criterion 12911 , Ac- 
cording to this criterion, in the case of 2 x 2 and 2x3 sys- 
tems, a density matrix gAB is inseparable (i.e., entangled) iff 
{qab) Ta (Ta- partial transposition with respect to system 
A) is not positive. Tables HI and HH show the results for anti- 
ferromagnetic case (J > 0). The parameters in the tables are 



TABLE I: Inseparability conditions of the output states in the three 
different scenarios of cloning. 
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TABLE II: Inseparability conditions of the output states in the three 
different scenarios of cloning, at low and high temperature limits. 
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FIG. 4: Entanglement phase diagrams of input and output 
states (achieved from three different schemes of entanglement 
cloning/broadcasting introduced in the text), when a = 1/ s/2. The 
regions labeled by 1 are the regions in which entanglement exists, 
whilst the regions labeled by 2 indicate no-entanglement regions. 
This figure shows that for 7 < 7 C (T > T c ), depending on the value 
of e, we may or may not have entanglement. -y c is a decreasing func- 
tion of M. In other words, the area of region 1 increases when M 
increases, as well. This may imply the advantage of the entanglement 
doner Mm over the other entanglement broadcasting schemes. 



as follows: 



y/(MI -l){M+l-AMe) 
' 4M(1 - e) ' 

>/(M + l)(3M- 1-AMe) 
" 4M(1 - e) ' 

y/3M 2 +2M-1 
4M ' 



, ,M + l + 2y/M 2 + M, 
7c = ln ( — 1 h 



(4) 



ei(2) = 



3M-1 
(M-lT4MJ)(l + cosh 7 ) 
2M[1 ± sinh 7 =F 25(1 + cosh7)] ' 
V(3M(1 - e ) - 1) (M(l - e) + 1) 
4M(1 - e) 



where 5 = ay/1 ~ a 2 , CI = < a < 1, and C2 = 
— 1 ^ a ^ 0. When 7 — > 00 and M = M;;;, since 
(3M - 1)/4M > (1 - M)/2M, there exists an overlap be- 
tween the e-inequalities in the third and fourth sub-rows of Ta- 
blellll In this case, one should notice that for (1 — M)/2M < 
e < (3M - I)/ AM, clones are entangled if \c? - 1/2| < af 
or I a 2 — 1/2 1 > af . This removes the ambiguity in such 
cases. 

Tables HI and HH imply that in most temperature regions, the 
inseparability inequalities are not symmetric with respect to 
a — > —a. In other words — unlike the case of e = — de- 
pending on the sign of a, the parameter regions over which 



the cloned pairs are entangled may be different. Another im- 
portant point (see the second row of Table |I| is the existence 
of a critical temperature T c (oc 1 /j c ) beyond which the cloned 
pairs for some a regions, \a 2 — l/2\ a c , for all es are not 
entangled. 

Overall, by taking into account the behaviors of the upper 
and lower bounds of the inseparability inequalities we can find 
that in some temperature regions, in Table U (Table HP, there 
exist intervals of a 2 (e) in which the cloned pairs are sepa- 
rable. The length of these intervals decreases when M in- 
creases (recall that Mm > Ma > M;). Furthermore, for a 
given a 2 (e) at intermediate (two limits of) temperatures, the 
range of e (a 2 ) in which the clones are entangled increases 
when M increases as well. Indeed, for some temperature re- 
gions, in Table U (Table Hill there exist some a 2 (e) in which 
clones for all e (a in CI or C2) and all three Ms are entan- 
gled — e.g., see first sub-row of Table U or second and fifth 
sub-rows of Table [TT] These facts together with the entan- 
glement phase diagrams in Fig. |H whose regions show ex- 
istence of entanglement or its non-existence for a = 1/ \2, 
indicate advantage of entanglement doner M^, over the other 
cloning schemes. That is, the optimal entanglement doner has 
an advantage over other mentioned schemes of entanglement 
broadcasting in the sense of robustness against thermal noise. 

Conclusion. — We have studied the role of thermal noise 
in some quantum cloning schemes through a simple model of 
temperature effect on spin states at the input of the cloning 
machines. The performance of the cloning machines depends 
generally on the values of the thermal perturbation coefficient, 
the orbit of the original state on the Bloch sphere, as well as on 
the temperature. In addition, three scenarios of entanglement 
cloning of thermally diluted two-qubit states have been inves- 
tigated. Our analysis shows that the clones generated from 
non-local transformations, in particular those out of the opti- 
mal entanglement doner, remain entangled for wider regions 
of parameters. I.e., the optimal entanglement doner shows 
a relatively larger region of entanglement in the parameter 
space. This can be considered as an advantage of optimal en- 
tanglement doner over the other scenarios in the sense of ro- 
bustness against thermal perturbations. This statement, how- 
ever, is subject to the thermalization model we have used; so 
for a general conclusion a more detailed study is still needed. 
Our results may be of importance in practical implementations 
of quantum cloning in systems in which thermal effects are 
unavoidable, e.g., nuclear spin systems lfl9l Eoll . Indeed, the 
large e regime of our approach — when r c is of the same or- 
der of magnitude as Tdiss. — has already been experimentally 
realized in a different context ltl9||. This can be considered as 
a non-economic cloning process 113011 . 
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